APS/123-QED 



Ionization of hydrogen atoms by electron impact at leV, 0.5eV and 0.3eV above 

threshold 



J. N. Das and S. Paul 

Department of Applied Mathematics, University College of Science, 92 A. 



P. C. Road, Calcutta 700 006, India. 



O 
O 

(N 

C 

(N 



o ■ 



(N : 
> ■ 

vo : 

o 

in 

o 

m 

O 

o 

>>: 
43 : 

Oh' 



X 



K. ChakrabartO 

Department of Mathematics, Scottish Church College, 
1 & 3 Urquhart Square, Calcutta 700 009, India 
(Dated: February 2, 2008) 

We present here triple differential cross sections for ionization of hydrogen atoms by electron 
impact at leV, 0.5eV and 0.3eV energy above threshold, calculated in the hyperspherical partial 
wave theory. The results are in very good agreement with the available semiclassical results of Deb 
and Crothers [T^ for these energies. With this, we are able to demonstrate that the hyperspherical 
partial wave theory yields good cross sections from 30 eV [l2^1 down to near threshold for equal 
energy sharing kinematics. 

PACS numbers: 34.80.Dp, 34.50.Fa 



I. INTRODUCTION 

Over the last couple of years considerable progress has 
been made in the study of ionization of hydrogen atoms 
by electron impact, apparently the simplest three body 
coulomb problem in quantum scattering theory, at low 
energies. Still full understanding of this problem, partic- 
ularly near threshold, has not been achieved. 

Ionization near threshold was studied by Pan and 
Starace where they reported a distorted wave cal- 
culations of (e, 2e) process with H, He and other rare 
gas targets at excess energies 4 eV (above the re- 
spective ionization thresholds) and below for equal en- 
ergy sharing kinematics. For atomic hydrogen they 
reported results at excess energies 2 eV and 0.5 eV 
above threshold in the coplanar constant 9 a b geome- 
try (in which the angle 6 a b between the emerging elec- 
trons remain fixed) with 6 a t, — tt. Jones, Madi- 
son and Srivastava |2j also reported a distorted wave 
(e, 2e) calculation with atomic hydrogen and helium tar- 
get for equal energy sharing kinematics but different ge- 
ometries. Their results were in good qualitative agree- 
ment with the experiments at 2 eV above threshold for 
atomic hydrogen for constant 6 a b geometries. However, 
for the other geometries presented, there were consider- 
able deviations from the experimental results. 

For equal-energy-sharing kinematics McCurdy and co- 
workers made a break-through calculation in their ex- 
terior complex scaling (ECS) approach 0- EH 0| ■ Their 
results for 30eV, 25eV, 19.6eV and 17.6eV agree excel- 
lently with the measured results of Roder et al H, Q, @ • 
However, below 2eV, results of ECS theory are not yet 
available. Another sophisticated approach is the conver- 
gent close-coupling (CCC) theory of Bray and associates 
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III Eil UK < which works beautifully for many atomic scat- 
tering problems and reproduces ionization cross section 
results very satisfactorily above 2eV excess energy. For 
2eV excess energy their results differ approximately by a 
factor of two from the absolute measured values |f|. Be- 
low 2eV excess energy CCC results are also not available. 
Recent calculations of Das and co-workers 0] for equal 
energy sharing kinematics in the hyperspherical partial 
wave (HPW) theory also reproduced the experimental 
data ja, |1] quite satisfactorily. 

So far the hyperspherical partial wave theory has not 
been tested below 2eV excess energy. Deb and Crothers 
[l3| have reported a semiclassical calculation that gives 
very good cross section results for low energies of 4eV and 
2eV above threshold and also for energies leV, 0.5eV and 
0.3eV above threshold. This calculation encouraged us 
to test whether the hyperspherical partial wave theory 
works for excess energy below 2eV. Here we made such 
a calculation for excess energy of leV, 0.5eV and 0.3eV 
above the ionization threshold. We found that the HPW 
theory gives cross section results in very good agreement 
with the semi classical calculation of Deb and Crothers 
fl3| for the above energies. One only needs to increase 
the asymptotic range parameter i?^ (defined later) to 
sufficiently large values of several thousands of a.u. It is 
interesting to note here that the hyperspherical 7?.-matrix 
with semiclassical o utg oing waves (HIZM- SOW) calcula- 
tion of Selles et al 14], for the double photo-ionization 
of the Helium atom, also requires Roo values of several 
thousands of a.u. for converged results. 



II. HYPERSPHERICAL PARTIAL WAVE 
THEORY 



The hyperspherical partial wave theory has been de- 
scribed in details in 0, and briefly in 0, . In 
this approach we determine scattering amplitude from 
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the T-matrix element given by 



Tfi = 



(1) 



where $i is the initial state wave function, Vi is the in- 
teraction potential in this channel and ' is the ex- 
act two-electron continuum wave function with incoming 
boundary conditions in presence of the nucleus, which 
is considered infinitely heavy and stationary at the ori- 
gin. The scattering state <I 'j ' is determined by expand- 
ing it in terms of symmetrized hyperspherical harmon- 
ics 0] which are functions of five angular variables. 
The hyperradius and the angular variables are defined by 
R = Vr{ 2 + r 2 2 , a = arctan(r 2 /ri), f\ — (Oi,(fii), f 2 = 
(0 2 ,4> 2 ) and U = (a,fi,f 2 ) and set P = y 'pi 2 + p 2 2 , 
a = arctan(p 2 /pi), Pi = (Q Pl ,<t> Pl ), Vi = (® P n<\> P i) and 
u)q = (ao,Pi,p 2 ), Pi, p 2 being momenta of the two out- 
going electrons of energies E\ and E 2 , and positions T\ 
and f 2 . 

Different sets of radial waves with definite fi = 
(L,S,tt), (where L is the total angular momentum, S 
is total spin and tt is the parity) satisfy different sets of 
coupled equations each of the form 
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These equations are truncated to N = N mx which is 
the number of channels retained in the calculation for 
each fixed fi. The N mx number of independent solutions 
of the truncated equations, need to be determined over 
the interval (0, oo). These equations may then be solved 
in different alternative approaches. One possibility is to 
partition this interval into three subintervals (0, A), (A, 
Roa) and (i?oo, oo), A being of the order of a few atomic 
units and Roc being a point in the far asymptotic do- 
main. The solution in (0, A) is then constructed as in 



the R-matrix ^(| method and then continued to Roc us- 
ing Taylor's expansion method 0,0. Beyond Roc the 
solutions are known from series expansions |l2j |. This 
approach, however, suffers from pseudo resonance prob- 
lems as pointed out in Ref. [18J and hence this is not 
the one adopted here. Other possibilities include solv- 
ing the set of equations as a two point boundary value 
problem (since the radial wave function is known at ori- 
gin, and at Roo from series expansions) as in the ECS 
method This again would require more computa- 
tional resources than that we have at present. The most 
effective approach for our purposes turns out to be the 
following. We construct N mx independent solutions of 
Eq. (2) over the interval (0, A) by solving these as a two- 
point boundary value problem. The k th solution vector 
is made to vanish at the origin and takes the k*' 1 column 
of the N mx x N mx unit matrix as its value at A. These 
solutions are then continued over (A, Roc) by solving 
the coupled system of equations by the Taylor's expan- 
sion method with frequent stabilization [2(j • Beyond Roc 
the solution may be obtained from expansion in inverse 
powers of R with suitable sine and cosine factors [T^ . [THI |. 
The asymptotic incoming boundary conditions then com- 
pletely define 0, the scattering-state wave function 
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y s . For the initial interval ( 0, A ) solution by the finite 
difference method proves most effective. In our earlier 
calculation a t higher energies we used a five-point 
difference scheme. This gives us very good cross sections 
for 30eV, 25eV, 19.6eV and 17.6eV for various kinematic 
conditions. In our present calculation we propose to use 
larger mesh size (double that of our previous calculation) 
and hence, in place of a five-point difference scheme we 
use a seven-point difference scheme. In the seven-point 
scheme we divide the interval [0, A] into m subintervals 
by using mesh points i?o, Ri, R 2l • • • , Rm-i, Rm where 
Rk = hk,(k = 0, 1,2, ■••m) and h — A/m. In this 
scheme we use the following formulae: 



J 



" 1 1 3 3 49 3 3 

fN( R k) = "72 [ — fN(Rk-3) ~ ^16/Af(i?fc-2) + -/jv(i?fc-l) - — /jv(-Rfc) + ^fN(Rk+l) ~ —l6f N (R k+2 ) 



(3) 



r 



for k = 3, 4, • • • , to— 4, m— 3. For k = 1, 2 and m— 2, to— 1 we use the the formulae 
I 
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h 2 
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M#o) + 6M#i) 



y h 2 f N (R 2 ) - ^f N (R 3 ) - h 2 f' N (R 3 ) + 6f N (R 4 ) + p N (RA) 
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(4) 
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f' N (R2) = i[ lf N (Ri) + 6fN(R 2 ) ~ y h 2 f' N (R 3 ) - ^f N (R 3 ) - h 2 f' N (R A ) + 6f N (R 4 ) + p N (R 5 ) } 

+ { —hS&^Uts)}. (5) 

1 10080 J vw/ v ; 

/jv(-^m-2) = To [ o/jv(-Rm-5) + 6/jy(-Rm-4) — h 2 f N (R m ^ i ) —f N (R m _ 3 ) — h 2 f N (R m - 3 ) 

fi 8 4 2 

+ 6f N (Rm-2) + ^fN(R m ~l)}+{- T ^h 6 f^\U)} (6) 
fN(Rm-l) = T/o [ o/jv(-Rm-4) + 6/Ar(i? m _ 3 ) — h 2 f N (R m ^ 3 ) —fN{R m -2) — h 2 f N (R m - 2 ) 

h z 8 4 2 

Q OQ 

+ 6/ Ar ( J R m _ 1 ) + -/ J v( J R m )] + {- IM ^/ l 6 / (Ml4) te)}. (7) 



In each of Eqs. (3)-(7) quantities on the right hand sides 
within the curly brackets represent the error terms. The 
corresponding difference equations are obtained by sub- 
stituting the values of second order derivatives from the 
differential equation (2) into these expressions. For con- 
tinuing these solutions in the domain (A, R^) we need 
first order derivatives f' N (R) at A. These are computed 
from the difference formula 



and hence we chose the above seven-point scheme in the 
present calculation. We included states with L up to 5. 
Values of L above 5 give insignificant contributions. The 
(Zi, I2) pairs which have been included in our calculation 
are sufficient for convergence as found from the results of 
calculations with the inclusion of larger number of chan- 
nels. All the computations were carried out on a 2 CPU 
SUN Enterprise 450 system with 512 MB RAM. 



- 128Mi? m -l) + 105/Ar(i? m )] 

+ yV^JU + i-^/^CO} (8) 

Here too, the quantity within curly brackets represents 
the error term. The solutions thus obtained in (0, A) 
are then continued over (A, i?oo) by Taylor's expansion 
method, as stated earlier, with stabilization after suitable 
steps 0| . 

III. RESULTS 

In our present calculation for the equal-energy-sharing 
kinematics and leV, 0.5eV and 0.3eV excess energies, we 
have included 30 channels and have chosen A = 5 a.u. 
(as in our previous calculation |12| for higher energies). 
Small variation of the value of A about the value chosen 
does not change the results. Here we need to choose Roo 
equal to 1000 a.u. for leV, 2000 a.u. for 0.5eV and 3000 
a.u. for 0.3eV for smooth convergence of the asymptotic 
series solution and for smooth fitting with the asymptotic 
solution |12| in the equal energy sharing cases. For un- 
equal energy sharing cases one may need to move to still 
larger distances. For going that far in the asymptotic 
domain a larger value of h (grid spacing) is preferable. 
Here we have chosen h = 0.1 a.u. up to A and a value 
0.2 a.u. beyond A in all the cases. Accordingly a seven 
point scheme, as described above, is more suitable over 
a five point scheme used in our earlier calculation 



For each incident energy three sets of results with dif- 
ferent geometry have been presented for the two outgoing 
electrons having equal energies. These are the constant 
- 9 a b geometry results, # a -constant geometry results and 
the results for symmetric geometry. For leV excess en- 
ergy we have presented these results in figures 1(a), 1(b) 
and 1(c) respectively. The corresponding results for en- 
ergy 0.5 eV are presented in figures 2(a), 2(b) and 2(c) 
respectively and the results for 0.3 eV have been pre- 
sented in figures 3(a), 3(b) and 3(c) respectively. The 
bottom row of each of the figures 1(a), 2(a) and 3(a) cor- 
responding to 9nh = 150° and 9 a b — 180° are as in our 
earlier work [21| (though these are now calculated with 
different R^ values). This is merely to ensure complete- 
ness in the results presented and to compare our results 
with the semiclassical calculation of Deb and Crothers. 
For other geometries, unfortunately, neither experimen- 
tal nor any theoretical results are available for compari- 
son. The overall agreement between our results and those 
of Deb and Crothers for 6 ab = 180° and 6 ab = 150° 
is highly encouraging. A little steeper rise of our results 
compared to those of Deb and Crothers ^| at 0° are in 
conformity with the general trends of our corresponding 
earlier results |l2| at 2eV and 4eV excess energies. Our 
results for 9^ = 120° and 100° also appear reasonable 
when compared with the shapes of the corresponding re- 
sults for 2 eV and 4 eV excess energy [l2j • Unfortunately 
there are no experimental results for verification. The 
results for (9 -constant geometry and those of symmetric 
geometry are also in very nice agreement in shapes, par- 
ticularly for 9 a = —30° and 9 a = —150°, with those for 
2eV and 4eV excess energy cases (see Das et al [T^l. 
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IV. CONCLUSION 

From the results presented above it appears that the 
hyperspherical partial wave theory works satisfactorily at 
leV, 0.5eV and 0.3eV excess energies. We have already 
good results 01 f° r energies up to 30eV for various kine- 
matic conditions. Calculations at a higher incident en- 
ergy of 54.4eV are now in progress and we propose to 
present them in future. Another point to note is that in 
this approach exploration of the far asymptotic domain is 



possible by increasing to thousands of atomic units. 
All these suggest that the hyperspherical partial wave 
theory is capable of being developed into a successful 
method for (e, 2e) collisions. 
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Figure Captions 

Figure 1(a). TDCS for coplanar equal-energy-sharing constant <d a b geometry at leV excess energy above threshold. 
Continuous curves, present results ; dashed-curves, semiclassical results of Deb and Crothers 

Figure 1(b). TDCS for coplanar equal-energy-sharing geometry at leV excess energy above threshold for fixed 8 a 
and variable 9b of the out going electrons. 

Figure 1(c). TDCS for coplanar equal-energy-sharing with two electrons emerging on opposite sides of the direction 
of the incident electron with equal angle 9 a at leV excess energy above threshold. 

Figure 2(a). TDCS for coplanar equal-energy-sharing constant O a (, geometry at 0.5eV excess e nerg y above 
threshold. Continuous curves, present results ; dashed-curves, semiclassical results of Deb and Crothers [Rj 

Figure 2(b). TDCS for coplanar equal-energy-sharing geometry at 0.5eV excess energy above threshold for fixed 9 a 
and variable 6^ of the out going electrons. 

Figure 2(c). TDCS for coplanar equal-energy-sharing with two electrons emerging on opposite sides of the direction 
of the incident electron with equal angle 9 a at 0.5eV excess energy above threshold. 

Figure 3(a). TDCS for coplanar equal-energy-sharing constant 8 a & geometry at 0.3eV excess e nerg y above 
threshold. Continuous curves, present results ; dashed-curves, semiclassical results of Deb and Crothers [Rj 
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Figure 3(b). TDCS for coplanar equal-energy-sharing geometry at 0.3eV excess energy above threshold for fixed 8 a 
and variable 6b of the out going electrons. 

Figure 3(c). TDCS for coplanar equal-energy-sharing with two electrons emerging on opposite sides of the direction 
of the incident electron with equal angle 6 a at 0.3eV excess energy above threshold. 
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Figure 2(c) 
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